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The recent experimental realization of ID equal Rashba-Dresselhaus spin-orbit coupling (ERD- 
SOC) for cold atoms provide a disorder-free platform for the implementation and observation of 
Majorana fermions (MFs), similar to the well studied solid state nanowire/superconductor het¬ 
erostructures. However, the corresponding ID chains of cold atoms possess strong quantum fluctua¬ 
tion, which may destroy the superfluids and MFs. In this Letter, we show that such ID topological 
chains with MFs may be on demand generated in a 2D non-topological optical lattice with ID ERD- 
SOC by modifying local potentials on target locations using experimentally already implemented 
atomic gas microscopes or patterned (e.g., double or triple well) optical lattices. All ingredients in 
our scheme have been experimentally realized and the combination of them may pave the way for 
the experimental observation of MFs in a clean system. 

PACS numbers: 03.75.Ss, 05.30.Fk, 03.65.Vf, 67.85.-d 


Majorana fermions (MFs) (lj obey non-Abelian ex¬ 
change statistics and are crucial for realizing fault- 
tolerant topological quantum computation HH3- Follow¬ 
ing initial theoretical proposals la? -Q, some possible 
signatures of MFs have been observed recently in ex¬ 
periments 15-21] using ID nanowires or ferromagnetic 
atomic chains on top of an s-wave superconductor and 
with strong spin-orbit coupling (SOC). However, these 
signatures are not conclusive because of disorder and 
other complications in solid state [2^ - l27j . In this con¬ 


text, the recent experimental realization of SOC L 28-33| 
in ultra-cold atomic gases provides a disorder-free and 
highly controllable platform for observing MFs. In exper¬ 
iments, ID equal Rashba-Dresselhaus SOC (ERD-SOC) 
and tunable Zeeman field have been achieved, which, to¬ 
gether with s-wave superfluidity, make it possible to ob¬ 
serve MFs 3414411] in ID atomic tubes or chains, similar 
as the nanowire systems. 


However, unlike solid state nanowire systems where s- 
wave superconducting pairs are induced from proximity 
effects, the superfluid pairing in the ID atomic chain is 
self-generated from the s-wave contact interaction, lead¬ 
ing to the strong quantum fluctuation that renders the 
long range superfluid order impossible in the thermody¬ 
namic limit. To circumvent this obstacle, quasi-ID sys¬ 
tems with multiple weakly coupled uniform chains [42l — 
[Hoj have been studied in both solid state and cold atoms, 
where transverse tunneling was found to lift the zero en¬ 
ergy degeneracy of multiple MFs. 


In this Letter, we consider a truly 2D non-topological 
fermionic optical lattice with the experimentally realized 
ID ERD-SOC. We raise the question whether single or 
multiple topological ID chains supporting MFs can be 
on demand generated at target locations in such non- 
topological 2D systems. Generally, a ID chain in a 2D 
lattice can be locally modified to satisfy the topological 
condition for MFs using the recently experimentally real¬ 
ized single site addressing (the atomic gas microscopes) 


5l|-55j or patterned (e.g., double or triple well) opti¬ 
cal lattices [56 h 58I ]. However, the atom chain is strongly 
coupled with neighboring chains through transverse tun¬ 
neling in the 2D system, therefore a naive expectation 
is that the coupling may destroy the local topological 
properties and the associated MFs. 

Here we show that ID topological chains with MFs can 
indeed be generated on demand from 2D non-topological 
fermionic optical lattices with experimentally achieved 
ID ERD-SOC. Local addressing lasers in atomic gas mi¬ 
croscopes can modify the effective local chemical poten¬ 
tials along single or multiple ID chains, leading to a 
quantum transition to discrete topological chains that 
are characterized by non-zero winding numbers and host 
MFs at chain ends. Multiple MFs in spatially separated 
multiple topological chains still couple, with the coupling 
induced zero energy splitting exponentially decaying with 
the distance of neighboring topological chains. Note that 
similar results apply also to 3D if ID chains can be locally 
addressed in a 3D optical lattice. In the weak trans¬ 
verse tunneling region (quasi-ID), the MF coupling is 
extremely small for two topological chains separated by 
one or two non-topological chains, making it possible to 
observe multiple MFs in 2D or 3D double or triple well 
optical lattices without requiring the single site address¬ 
ing. 

Model: We consider a spin-1/2 ultra-cold degenerate 
fermionic gas (spin j" and j.) in a 2D square lattice with 
the lattice size N = n x x n y . As shown in the schematic 
picture Fig. H] two Raman lasers couple two spin states to 
induce ID ERD-SOC along the x-axis. The far-detuned 
local addressing lasers [5ll455| can modify the local po¬ 
tential of the optical lattice along a ID chain at target 
locations along the ^-direction. Multiple spatially well 
separated ID chains can also be generated using addi¬ 
tional local addressing lasers. In the 2D lattice, the tight- 
binding mean-field BdG Hamiltonian 


Hb<ig = H l + Hsoc + Ho + Ha, 


(1) 
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FIG. 1: Illustration of the proposed experimental setup. Grey 
arrows represents 2D square optical lattice lasers. Red tube 
demonstrates the ID potential chain induced by local address¬ 
ing lasers in atomic microscopes (red arrows). Two counter- 
propagating Raman lasers (oran ge arrows) couple two spin 
states, generating ID ERD-SOC I28l - f33ll . 


where Hl 

Ei^c^ci. 


— Xa.o-,7;^ (^i }lT Ci+r/,cr + h.c.'j 

is the bare Hamiltonian in the 2D 


lattice with rj = {x,y}. The fermionic operator Cj 
(Ci iCT ) creates (annihilates) a particle with spin a at 
site i = (i x , i y ). We use p = p — 2 t x — 2 t y for the 
effective chemical potential to match with that in the 
continuous model. t x and t y are the nearest neighbor 
tunnelings along x and transverse y directions respec¬ 
tively. H S oc = “ C^Ci+xd + H - c -) + 

h z X)i(Ci -fCi.t — Clfk.i) describes the experimentally 
available ID ERD-SOC, with the SOC strength a and 
Zeeman field h z . Hjj = JA crVr{iy)Cj a Ci }Cr represents 
the ID potential dip with the local potential Vr{i y ), 
which is generated by local addressing lasers and non¬ 
zero only at the i y chains. H& = — JA Ai (Cj^Cj^ + h.c.) 
is the mean-field paring Hamiltonian, with the order 
parameter Ai = —g(CipCi^) and on-site interaction 
strength g. Hereafter we take t x = t as the energy 
unit. We solve the corresponding BdG equation self- 
consistently with the pairing gap equation and fixed 
chemical potential, following the standard numerical 
procedure |H-[63l. The order parameter is chosen to 
be complex to find the ground state. We use the box 
boundary condition for the self-consistent calculation. 
In practical experiment, there is a weak harmonic 
confinement which may change the (location of MFs, but 
don’t change the essential physics 63]. 

The above BdG Hamiltonian preserves particle-hole 
symmetry 5IT BdG S _1 = - H Bc lg , where 5 = t x IC, 
t x = t x <g> do (B> po, Ti, di are 2x2 Pauli matrices acting on 
particle-hole and spin spaces respectively, po is a N x N 
identity matrix on the lattice site space, and K, is the 
complex conjugate operator. If the order parameter Aj is 
real, the Hamiltonian is also real, which preserves a time- 



FIG. 2: MFs in single ID topological chain in the 2D optical 
lattice, (a) The amplitude of the order parameter |A|. (b) 
The quasiparticle energy spectrum, (c) The zero energy mode 
wave function along the central chain, (d) Zero energy mode 
density (Log scale). Parameters: t v = t, a = 2t, q = —5.5 1, 
h z = lAt, V = -1.45f, p = -1.555 1. 


reversal like symmetry QHbcLG® 1 = HbiIG with 0 = /C, 
as well as a chiral symmetry SHBdGS _1 = —Hb<ig with 
S = 0 • 2 = t x . In this case, the system belongs to the 
BDI topological class characterized by a Z topological 


invariant 


M 


One topological chain: First consider a 2D lattice with 
no tunneling along the y -axis (t y = 0), thus the 2D 
lattice is composed of individual x-direction ID chains. 
At the central chain we add an extra constant potential 
V T (y c ) = V, so that the central chain becomes topologi¬ 
cal, while other parts of the system are still in the non- 
topological region. Here the topological region is defined 
locally by the criteria h z ^ \/(p — Vt) 2 + A 2 , same as 
the usual ID topological chains [llj,[l4j. In our numerical 
calculations, we take n x = 81, n y = 9, and y c = 5. In 
the central topological chain, two MFs should exist at two 
ends. When t y is increased from zero to t, the topologi¬ 
cal chain couples with neighboring non-topological chains 
and the system changes from ID to quasi-ID and finally 
to 2D. A nature question is whether MFs at the center 
chain will survive with the strong coupling. 

Fig. [2] demonstrates the existence of MFs even in truly 
2D region with t v = t. The amplitude of the superfluid 
order parameter Ai is plotted in Fig. [2]]a). We find that 
Ai is homogeneous along the z-axis in the self-consistent 
calculation, except at the boundary. Furthermore, Ai has 
a constant phase across on the whole 2D system, there¬ 
fore we can choose it to be real without loss of general¬ 
ity for the discussion of the topological properties. Fig. 
(2l]b) shows the quasiparticle energy spectrum, where we 
clearly see the existence of Majorana zero energy modes 
with a tiny energy splitting E « ±2 x 10 _5 t mainly due 
to the finite size effect. The mini-gap energy, defined as 
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the energy difference between the zero energy mode and 
the next lowest quasiparticle state, is comparable to the 
order parameter Ai. Fig. [^c) shows the wavefunction 
of the zero energy mode (E ss +2 x 10 _5 f state) along 
the central chain, which satisfies the relation for MFs: 
u a = Xv ai X = ±1, indicating the central chain is still 
topological with two MFs at its ends. Fig. [2] (d) shows 
the density of the zero energy mode, which is square of 
the zero energy mode wave function for both spin up 
and spin down atoms in the 2D plane. The zero energy 
mode still localizes at the ends of the central chain, but 
slightly spread to neighboring chains which is due to the 
transverse tunneling. Note that in practical experiments, 
there exists a finite detuning for the Raman coupling be¬ 
tween two bare spin states, which corresponds to an in¬ 
plane Zeeman field h y a y in our notation. Such non-zero 
in-plane Zeeman field is known to break the inversion 
symmetry and lead to the FF type of ground states with 
finite momentum pairing. We confirm that our results 
still hold in the FF state (see Supplemental Material). 

Topological characterization: The emergence of MFs 
at the edges of the central chain originates from the 
bulk topological properties of the 2D optical lattice with 
the imprinted ID topological chain. In the above self- 
consistent BdG calculations, both order parameter and 
atom density are almost homogeneous along the z-axis, 
therefore it would be a good approximation to assume 
that the bulk is uniform. With a periodic boundary con¬ 
dition along the x-axis, the momentum k x is a good quan¬ 
tum number. The 2D lattices can be taken as a series of 
individual ID chains coupled through transverse tunnel¬ 
ing, with the effective BdG Hamiltonian 


H B dG{k x ) = H 0 (k x )p 0 + {Vt z ct o + A'T y cr y )p' - t y T z a 0 p x , 

( 2 ) 

where Ho(k x ) = [—2t x cosk x — p]t z oq + 2asin k x r z a y + 
h z T z a z +AoT y a y describes the original uniform individual 
chains with the SOC and Zeeman field, p spans the y- 
axis chain space with po as the identity matrix, (p')ij = 1 
for i = j = y c and 0 otherwise. The p' part describes the 
potential and order parameter differences of the central 
chain from others. The p x term describes the i/-axis hop¬ 
ping between nearest neighboring chains, with (p x )ij = 1 
for |z — j | = 1 and 0 for others (i,j= 1 • • • n y ). 

The topological properties of the BdG Hamiltonian ([2]) 
can be characterized by the winding number W. For 
a single ID chain in the 2D optical lattice, the BdG 
Hamiltonian is in the BDI topological class with a chiral 
symmetry SHBddk^S -1 = -H BdG (k x ), where S = t x . 
Therefore the BdG Hamiltonian can be transformed to 
be off-diagonal in the r space 


UH BdG U+ 


0 A{k x ) 
A T (-k x ) 0 


= h{k x )r x + d,T y (3) 


by a unitary transformation U = e i G/ 4 )' r H j where 
A(k x ) = h{k x ) - id , h(k x ) = {[-2t x cosk x - p]a 0 + 



FIG. 3: Complex value z as function of k x in the region (0, n) 
with (a) V = — 0.8t (non-topological) and (b) V = —1.2 1 
(topological). The other parameters: t y = t, a = 2 1, Ao = 
0.4f, h z = lAt, A' = O.lf, p = —1.555 1. (c) Band structure 
at the topological phase transition point V = —t. 


2asmk x a y +h z a z }p 0 -t y o 0 p x +Vo 0 p ', and d = A 0 a y p 0 + 
A 'a y p'. The winding number 


W = - 


i 


7r 



(4) 


where z(k x ) = det A{k x )/\ det A(k x )\ 166l ]. 

Fig. [3] shows the change of the topological properties 
with the potential V along the central chain. When V is 
small, the whole lattice, including the center chain, is in 
the non-topological region. The corresponding complex 
value of z rotates when k x changes from 0 to k x = 7r 
as shown in Fig. [3] (a), indicating W = 0. When the 
potential depth \V\ increases beyond a threshold value 
V c ks — t, W changes to —1. Across V c , the band gap 
closes (Fig. [3] (c)) and reopens, indicating a topological 
phase transition to a phase where the central chain be¬ 
comes topological and hosts a pair of Majorana fermions, 
agreeing with the self-consistent calculation. 

Multiple topological chains: Multiple topological 
chains may be generated using additional local address¬ 
ing lasers to obtain multiple MFs. We first consider 
two topological chains separated by one non-topological 
chain, with the extra potential Vr adding at y = 4 and 
y = 6. The energy spectrum from the self-consistent 
BdG calculation is plotted in Fig. |4](a), showing one zero 
energy mode although there are two topological chains. 
This is due to the strong coupling between two chains, 
leading to the winding number W = — 1, instead of —2. 
Therefore there is only one pair of MFs when two chains 
are close. Fig. [4] (b) shows the density of the zero energy 
mode, which widely spreads along the y- axis from y = 4 
to y = 6. 

When two topological chains are further separated by 
more than one non-topological chains, the analysis of the 
bulk topological properties based on constant order pa¬ 
rameter phase shows that W = — 2, indicating two pairs 
of MFs. However, in the self-consistent calculation, the 
order parameter phase is no longer uniform due to the 
interaction between two MFs at the same end, leading to 
the splitting of the zero energy states. In Fig. [2(c), we 
plot the quasiparticle energy spectrum for two topologi- 
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FIG. 4: MFs in two topological chains at y = 4 and y = 6 
(a),(b) or at y = 2 and y = 8 (c),(d). (a) Quasiparticle energy 
spectrum for two topological chains separated by one non- 
topological chain, (b) Density of the zero energy mode (Log 
scale), (c) Quasiparticle energy spectrum for two topological 
chains separated by five non-topological chains. The inset 
plots the lowest positive quasiparticle energy as a function of 
the number of non-topological chains in between, (d) The 
phase of the order parameter A. t y = t, the other parameters 
are the same as Fig. [2] 


cal chains located at y = 2 and y = 8 obtained from the 
self-consistent calculation. Fig. []] (d)) shows the phase 
9{x,y ) of order parameter A = \/S\e l6 ^ x ' v \ which has an 
antisymmetric structure between two topological chains. 
The phase difference between two ends of one topological 
chain is opposite from the other topological chain. 

In principle there should be no zero energy modes left 
due to the interaction between MFs, which splits the en¬ 
ergy away from zero to a finite value 45, [46 1. In prac¬ 
tice, due to the large distance between two topological 
chains, the coupling strength between two MFs on dif¬ 
ferent topological chains is extremely small and the en¬ 
ergy splitting becomes negligible. For instance, the five 
chain separation in Fig. [4] (c) leads to an energy splitting 
E « 8 x 10 ~H. The inset in Fig. [4] (c) shows the change 
of the splitting with the number of non-topological chains 
between two topological ones. When the number is 0 and 
1, W = —1, and the splitting is almost zero since there 
is only one MF at each end. When the number is 2 and 
above, the winding number becomes —2 and the interac¬ 
tion of two MFs induces a splitting. The energy splitting 
decreases exponentially with the distance between two 
topological chains and approaches almost zero for the 5 
lattice separation. 

Multiple MFs in superlattices: The interaction between 
MFs in topological chains can also be significantly re¬ 
duced by decreasing the tunneling along the transverse 
direction, which makes the system quasi-ID, instead of 
2D. In this case, no large separation between neighboring 
topological chains is needed, making it possible to gen¬ 


FIG. 5: MFs in quasi-ID triple-well superlattices, (a) Illustra¬ 
tion of the triple-well superlattices along the j/-axis. Topolog¬ 
ical chains (TC) are at y = 2, y = 5 and y = 8, separated by 
non-topological chains (NTC). (b) Quasiparticle energy spec¬ 
trum showing three pairs of zero energy modes, (c) Density of 
zero energy modes (Log scale). Other parameters t y = O.lt, 
a = 0.75f, g = -3.5f, h z = 0.7f, V = -0.55 1, p = -0.55 1. 


erate multiple MFs using patterned optical superlattices 
along the y-axis. In experiments, optical superlattices 
such as double well or triple well lattices can be generated 
using the superposition of different lattice beams (fhi - IHsl , 
which are much easier than the single site addressing. 
Fig. [5] (a) shows a triple well optical lattice along the 
y -axis with one of the triple wells in the topological re¬ 
gion (i.e., two neighboring topological chains separated 
by two non-topological chains). With a small transverse 
tunneling t y = O.lt, the energy splitting for the MF zero 
energy state is as tiny as E ss 5 x 10~ 5 f, as shown in Fig. 
[5] In our calculation, we put three topological chains 
at y = 2, y = 5 and y = 8, and find one pairs of MFs 
formed at each topological chain ends. We also confirm 
that similar physics occurs if another triple well superlat¬ 
tice is applied along the 2 -axis to form a 3D lattice with 
weak tunneling along both y and 2 directions. 

Conclusion: We show that, with the assistant of 
atomic gas microscopes or patterned optical super¬ 
lattices, 2D non-topological optical lattices with exper¬ 
imentally achieved ID ERD-SOC can host non-coupled 
ID topological chains with MFs. We emphasize that al¬ 
though we illustrate the idea using a 2D geometry, the 
same physics also applies to 3D non-topological optical 
lattices, providing selected ID chains can be locally mod¬ 
ified or patterned optical superlattices are used. All in¬ 
gredients in our proposed schemes are available in current 
experiments, and the scheme may lead to unambitious ex¬ 
perimental signature of the long-sought MFs in a clean 
cold atomic system. 
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Supplemental Material 
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FIG. 6: (a) Density of the zero energy mode (Log scale), 
(b). Quasiparticle energy spectrum. h\\ = 0.2£, the other 
parameters are the same as Fig. 2 in the main part of the 
paper. 


In this Supplemental Material, we show the existence 
of MFs in a single ID topological chain in the 2D op¬ 
tical lattice with an additional in-plane Zeeman field 
h\\. HsoC = ~ + H.c.) + 

h z £i(C'? t Ci,t-C^CU) + /dl £i (iC+ t C iA + H.c.), where 
the last term represents the in-plane Zeeman field. The 
self-consistent BdG results are present in Fig. [G] 








